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Abstract 



We show that a cosmic string associated with spontaneous U{1)r symmetry breaking gives 
a constraint for supersymmetric model building. In some models, the string can be viewed 
as a tube-like domain wall with a winding number interpolating a false vacuum and a true 
vacuum. Such string causes inhomogeneous decay of the false vacuum to the true vacuum 
via rapid expansion of the radius of the tube and hence its formation would be inconsistent 
with the present Universe. However, we demonstrate that there exist metastable solutions 
which do not expand rapidly. Furthermore, when the true vacua are degenerate, the structure 
inside the tube becomes involved. As an example, we show a "bamboo"-like solution, which 
suggests a possibility observing an information of true vacua from outside of the tube through 
the shape and the tension of the tube. 



1 Introduction 



The global U(1)r symmetry plays an important role in supersymmetric field theories, in par- 
ticular in supersymmetry (SUSY) breaking [IHS] (See [6H8] for reviews and references therein). 
In [2l[3l[9], by exploiting the Nelson-Seiberg theorem [I], a connection between metastability 
and R-symmetry was demonstrated in the context of generalized Wess-Zumino models with 
generic superpotential. From more phenomenological viewpoint, the U(1)r symmetry must 
be broken explicitly or spontaneously to generate Majorana gaugino masses. Gaugino masses 
are not induced by SUSY breaking without the U(1)r symmetry breaking. 

Indeed, several types of models for the U(1)r symmetry breaking have been studied [TU1 - 
[23]. In some models, the vacuum with both SUSY and U(1)r breaking may be a global 
minimum. However, in many models such vacuum is a metastable minimum and there is a 
global minimum, where SUSY and U(1)r may be unbroken. 

Through the cosmological phase transition, there may appear solitonic objects such as 
domain walls, cosmic strings and monopoles [25] through the Kibble-Zurek mechanism [261127] . 
When a global U(l) symmetry is spontaneously broken, a global string appears [28]. Thus, 
when the U(1)r symmetry is broken spontaneously in SUSY models, there would appear a 
global string, which we refer as an R-stringJ^ The R-string would be stable in those models 
in which the U(1)r breaking vacuum is a global minimum, and that would lead to several 
cosmologically interesting aspects [3Tj . 

On the other hand, when the U(1)r breaking vacuum is metastable and the model has 
another global minimum with SUSY and U(1)r unbroken, there may appear an R-string, 
whose core corresponds to the true SUSY vacuum, i.e., R-tube. One may think that such 
an R-tube is unstable because the energy density in the core, which is the SUSY vacuum, 
is lower than one outside, which is the SUSY breaking metastable vacuum. Thus, it would 
"roll-over" and the true SUSY vacuum would expand in the Universe [32HM] . In this case, 
the SUSY breaking could not be realized successfully. One may conclude that a scenario with 
R-tube formation is ruled out by this mechanism. 

However, since the domain wall tension works as a centripetal force for R-tube, its radius 
may be stabilized if the domain wall tension is large enough and the energy discrepancy 
between SUSY vacuum and SUSY-breaking vacuum is small enough. In such a case, the 
cosmological disaster can be avoided. The (in) stability of the R-tube soliton depends on 
parameters in the SUSY models. In principle, we can have constraints on SUSY-breaking 
models from this consideraion because (in) stability of the R-tube is determined by parameters 

1 See also for another type of strings, which appear through SUSY breaking [291130] . 
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of SUSY breaking models. Note that such constraints are independent of the requirement 
that the metastable vacuum decays slowly into the true vacuum by the tunneling effect [35] , 
compared with the Universe age. Therefore, it is quite important to study the R-string/R-tube 
formation and its (in) stability. Some relevant studies have been carried out in Refs. [3S1EZ]- 

In this paper, we study in detail the structure of the R-string/tube solution in SUSY 
models. In a simple but (semi)realistic SUSY breaking model, we study stability of the 
R-tube by exploiting a piecewise linear approximation and numerical solutions. By using 
linear approximation, we obtain constraints for the stable R-tube. We also show examples 
of (meta)stable/unstable R-tube configurations numerically. We emphasize that the winding 
number, which is an important quantity to characterize features of the R-tube solutions, is 
also relevant to the stability of the R-tube. 

We also show that the core of the R-tube can have more complicated structure in certain 
SUSY models where the true SUSY breaking vacua are degenerate. For example, suppose 
that the SUSY model has a Z 2 symmetry and it is broken at the true SUSY vacuum. Then, 
the core of the R-tube would be separated into two vacua by a domain wall. Since it looks 
like a (gourd-shaped) bamboo, we refer it as the bamboo solution. We also study aspects of 
the bamboo solution. Other types of structure inside the core of strings would be possible. 

This paper is organized as follows. In section 2, we illustrate the R-string and tube 
solutions in simple models as warm-up. In section 3, we study the R-tube in a (semi)realistic 
but simple SUSY breaking model, that is, an O'Raifeartaigh-like model with non-canonical 
Kahler metric. We analyze the (in) stability of the R-tube numerically at several parameters. 
In section 4, by showing the bamboo solution, we demonstrate the fact that quantum number 
in the SUSY vacua significantly affects the shape and the tension of the string. Section 5 
is devoted to conclusion and discussion. In appendix A, we show basics of the relaxation 
method for solving a differential equation. 

2 Stable R-string and tube solutions 

Before going to detailed studies of metastable strings, which will be shown in the next section, 
we would like to present simple stable solutions as a warm-up. Here, we illustrate a stable 
R-string and a stable R-tube which is a tube-like domain wall with winding number, by using 
single complex scalar field models. 
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2.1 R-string 

Consider the following simplest spontaneous R-symmetry breaking model. Superpotential is 
linear in a chiral superfield X which will be a trigger for SUSY breaking, 

W = fX. 

To stabilize the pseudo-moduli X in this SUSY breaking vacuum, we introduce the following 
non-canonical Kahler potential by hand, 
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xx 



i- T \x\ 



(2.1) 



Thus, the potential of this theory is given by 



v(x) = f-^ x \xn^\x\\ 



(2.2) 



This model can be viewed as a low energy effective theory of one of the O'Raifeartaigh 
models studied in [12]: When the pseudo- moduli space is stable everywhere along messenger 
directions, by integrating out the messengers, one obtains non-trivial corrections to the Kahler 
potential. Expanding the Kahler potential up to (9(|X| 6 ), one can reproduce a theory similar 
to ( 12. 21) . On the other hand, when the pseudo-moduli space has a tachyonic direction at 
a point in the space, which is phenomenologically interesting situation in gauge mediation 
models [3H2D], the existence of messengers is crucial and two-field model is required. This is 
the main topic in the next section. 

When n 2 x > and Ax > 0, the X field develops its vacuum expectation value and the R- 
symmetry is broken (X has the R charge 2). The minimum of the potential V(X) is obtained 
at \X\ = X m - m = a/2/U^/Ax- Note that this vacuum is the global minimum of the potential 
V(X). Since the global U(1)r symmetry is spontaneously broken, the R-string would be 
formed. 

Let us introduce dimensionless variables as 



X X m - m T, x. 



2/4 s 



Ax/ 



2 X ^ 



Then the effective Lagrangian is given by 



A = ^>0. 
4/4 



(2.3) 



V(T) 



, V(T) = 1 
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(2.4) 
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Figure 1: The minimal winding solution of the R-string (solid line) for A = 1/2: the profile 
function is shown in the left panel and the corresponding energy density is shown in the right 
panel. The broken lines stand for the standard global vortex which is the solution of the 
model with the minimal kinetic term in Eq. I\2.4\) . 



Here a positive definite metric at the minimum (T = 1) requires A > 1/4. We use the 
dimensionless cylindrical coordinate (p, 9, z) for constructing a straight R-string along the 
z-axis. We make the following standard Ansatz, 

T(p,9,z) = f(p)e ine , (2.5) 

where f(0) = and f(p) — > 1 at p — > oo. We numerically solve the equation of motion for a 
minimal winding solution (n = 1). The solution is shown in Fig. [TJ 

For later convenience, let us estimate a size of R-string, R, by using the following simple 
approximation 

f{p) = (£ftorp<R and f(p) = 1 for p > R, (2.6) 

where the power of p is determined by requiring smoothness of the configuration at p = 0. 
The total energy of this configuration, E, per the string length Az is estimated as 



c nW + T77TT lQ g d + o V(1)A 2 + --^i? 2 , a n = — — — , (2.7) 



2//^2ttAz v ' V(l) to V^/ 2 24A ' l + 3n + 2n 2 

where an i? independent constant c n (A) which should be numerically determined is intro- 
duced, and an IR-cutoff A is also introduced in order to regularize a well-known logarithmic 
divergence of a global vortex. The above energy takes the minimum at R ~ Rstrmsi 

2n 2n 2 _ 1 / 1 \ 

-Rstring = — p= ~ , % = - I 1 - — , (2.8) 

Ja^m T m T A \ 4A/ 
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Figure 2: The transverse sizes of the R-strings, Rt, with the linear approximation (solid line) 
and numerical calculations (points) using the relaxation method reviewed in Appendix lAl 



which is the transverse size of the R-string. Here is the dimensionless mass of T in the 
vacuum T = 1. For instance for A = 1/2 and n = 1, we obtain ttit = 1 and -R s trmg = 2\/3. 

To check the approximation (12.61) by comparing with numerical calculations, we introduce 
another definition of a transverse size of the R-string as 



Rt = 



f™dpp 2 K T 
J™dppK T ' 



Kn 



(f(p)) 2 

v(f( P )y 



(2.9) 



Here Kt gives a finite contribution from the kinetic term along the p direction into the 
energy and is useful to define the transverse size. Note that this quantity does not include 
the cut-off dependence. We compute this both analytically with the linear approximation 
and numerically, see FigfSJ For instance we observe Rt = 1.95 for A = 1/2 by a numerical 
calculation. As can be seen in Fig. (2J the linear approximation nicely reproduces the numerical 
results (we need to pay attention to errors of 10% — 30% in the linear approximation). 

2.2 Tube solution 

Here, in order to illustrate the tube solution, we study a non- super symmetric bosonic theory 
as a toy model. Let us study the model with the following scalar potential, 

\X\ 2 



V(X) 



M 2 



2\2 



(2.10) 



and the canonical kinetic term, where M and v are taken to be real. This model has a 
global U(l) symmetry (no longer U(1)r symmetry), under which the field X transforms. 
This potential has two degenerate vacua, that is, |X| =0 and v. At the former vacuum, the 
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Figure 3: The minimal winding solution of the string of the model given in Eq. A2.12\) . The 
profile function is shown in the left panel and a slice at a hxed z of the corresponding energy 
density is shown in the right panel. 



global U(l) symmetry is unbroken, while the U(l) symmetry is broken at the latter vacuum. 
Then, a global string would be formed. Again, let us rescale the fields and coordinates as 



X = vT, Xf, = — 



then the Lagrangian becomes 



C 



M 2 



d^T 



2\2 



|T| 2 (l - \T\ 2 ) 



We make the Ansatz for the minimally winding string, 

T = f{r)e in \ 



(2.11; 



(2.12) 



(2.13) 



where f(0) = and f(r) — > 1 at r — > oo. The solution is again obtained numerically which is 
shown in Fig. [3J As can be seen in Fig. [31 the string has a substructure that is a hole inside 
the string. Thus, we refer it as the tube. It is the asymmetric phase (X ^ 0) outside the tube 
while it is the symmetric (X = 0) phase inside it. 

This tubelike string solution can be regarded as a tube of a domain wall. Indeed, there 
also exists a domain wall in this model. For instance, a solution interpolating the two vacua 
T = 1 at X\ = oo and T = at xi = — oo is given by 



T 



1 



VT 



-2xi 



(2.14) 
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Figure 4: Numerical results of a ratio Rx/Zn 2 for higher winding solutions, by using the 
relaxation method with efficiently large relaxation time, r > 20 x n 4 . 



with a dimensionless tension T wa u = 1/2. Thus, assuming that the field configuration of the 
tube along the radial direction is well described by this solution, the total energy E per the 
tube length Az of the tubelike solution with a radius p = R can be estimated by 

^5§7 ssT ">> fl + " 2|og (£)- (2 - 15 » 

as long as the "thickness" of the wall is much smaller than the radius R. Minimizing this, we 
get the transverse size of the tube solution 

n 2 

R « - — = 2n 2 . (2.16) 
T wa \i 

Note that the stabilization mechanism of the tube solution is different from that of the R- 
string (without a hole) where the kinetic energy and the potential energy are balanced. As a 
result, the transverse sizes have different dependences on the winding number n as R oc n for 
the R-strings and R oc n 2 for the tubes, respectively. 

We can define a transverse size Rt of this tube-like string similar to Eq. (12.91) with Kt = 
{f'{p)) 2 and observe Rt = 2.06 for the minimal winding tube by a numerical calculation. 
Ratios Rt/2u 2 with higher winding solutions are listed in Fig]U It suggests that the above 
approximation works well. 

3 R-tube and Vacuum Instability 

3.1 Metastable R-tube 

In section 12.11 we studied the single field SUSY breaking model as a toy model of one of 
the O'Raifeartaigh models discussed in [12] in which classical pseudo-moduli space is stable 
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everywhere. In this section, we move on to phenomenologically interesting situation where 
pseudo-moduli space has a tachyonic direction, in which large gaugino masses are generated 
by gauge mediation [2"U|I3"8]. In such models, the R-symmetry breaking vacuum is metastable, 
thus the R-string solution can be a tube-like domain wall with winding number as showed in 
section 12.21 

Here, we study the illustrating supersymmetric model with two superfields, X and <ft. 
These superfields have R-charges, R[X] = 2 and R[<f>] = 0, and the superpotential is given by 



W = X(f) 2 - /j 2 X. 



In addition, we consider the following effective Kahler metric, 



3xx 



1m? 



|X| 2 + 



A 



4m 4 



\X 
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X(j> 



hk = o. 



(3.1) 



(3.2) 



This model has the Z2 symmetry, under which X and <p are Z2 even and odd, respectively. 
This model has the discrete SUSY vacua, 



X = 0, 



(3.3) 



and the SUSY breaking vacuum, 



X 



m 

7x' 



(3.4) 



where the U(1)r symmetry is also broken. The former is the true vacuum, while the latter is 
the metastable vacuum whose vacuum energy is V = /i 4 V(l) = /i 4 (l — jr). 
For later convenience, let us introduce dimensionless variables by 



m 

then the Lagrangian is of the form 

1 



m 



■x 



m 



C 



V(T) 



V{T)\s< 



\TV\s\ 



(3.5) 



(3.6) 



This Lagrangian is characterized by two dimensionless parameters A and e. For instance, in 
the SUSY breaking vacuum (T, s) = (1, 0), dimensionless masses for T and s are wriiten by 



1 



ran 



A 



1 



T- 1--T- 



4A 



respectively, and that is, existence of the SUSY breaking vacuum requires 



1> Tx >0 - 



4A 



> 1. 



(3.7) 



(3.8) 
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Figure 5: The minimal winding solution of the R-tube of the two-scalar model for e = 1, 
A = 0.27 given in Eq. A3.6}) . The profile functions (solid line for \T\ and broken line for s) are 
shown in the left panel and a slice at a fixed z of the corresponding energy density is shown 
in the right panel. Here we chose the plus sign (7 = +1). 



If one is interested in vacuum selection, a simple criterion is a ratio of tachyonic masses 
at the origin (T, s) = (0, 0) where we have 

K),t) 2 = ~, (m ,,) 2 = -|. (3.9) 

Since in the early universe field values are assumed to be around the origin, if tachyonic 
mass of T is larger, one may expect that supersymmetry breaking model is preferable^. An 
inequality 

I < h < 3 - 10 > 

is required for selecting the SUSY breaking vacuum. 

Now we are ready to construct the R-tube in this two-scalar model. To this end, we make 
the Ansatz 

T = f(p)e me , s = 1 h(p), (3.11) 
with 7 = ±1. Because of the Z 2 symmetry, the solutions of f(p), h(p) are independent of 7. 



Similar to the model in section I2.2[ it is the symmetric phase inside the tube while it is the 
asymmetric phase outside the tube. A sharp contrast among two models is that the outside 



2 Here we have been studying a simple toy model. To discuss vacuum selection more seriously, one need to 
go back to an original realistic model and specify a history of the early universe. 
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is the true vacuum in section 12.21 and is the false vacuum in this section. One may guess that 
stable solution does not exist since the core of the tube has lower energy than its outside 
and hence larger radius would be favored energetically, which causes the "roll-over" problem. 
However, because the tension of the wall acts as a centripetal force for the R-tube, we will 
find that there exist metastable tube-like field configurations. In order to see a typical R-tube 
numerical solution in this model, here we show an example in Fig. [5] with e = 1, A = 0.27. 
Note that the profile function of the winding field T, whose mass is very small, has a very long 
tail compared to that of the solution in section 12.21 On the other hand, the unwinding field 
s, whose dimensionless mass is of order 1, converges exponentially. In the next subsection, 
we will investigate stability of the R-tube by varying those parameters. 



3.2 Instability of R-string and Broken Z 2 Symmetry 

If we set s = to keep the Z 2 symmetry, the model discussed in this section reduces to 
just the model discussed in section 12.11 except for an overall factor. Therefore the R-string 
solution (s = 0) without a hole inside is also a solution in this model. However, such an 
R-string would be almost always unstable and transforms into an R-tube with non-zero s 
inside. Since non-vanishing s means the broken Z 2 symmetry, we observe below that the Z 2 
symmetry inside the R-tube in this model is almost always broken. 



login!-*? 2 ] 



log l0 [-<r] 
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Figure 6: Tachyonic mass of s around R-string for A = 0.27 (the left panel) and A = l/2(the 
right panel). 



Let us consider an infinitesimal fluctuation 1) around the R-string solution discussed 
in section [2~T1 and study whether a direction along s is tachyonic or not. A linearized equation 
for s is given with an eigenvalue q 2 as 



-- {PS > 



--V(T) + -\T\ 



T=T a , 



q 2 s. 



(3.12) 
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For instance we observe tachyonic masses of s numerically for many sets of parameters {A, e} 
as shown in Figj6j We therefore make a conjecture that 

q 2 < 0, (3.13) 

for almost all the winding number n and the almost whole parameter region of {A, e} satisfying 
the inequalities (13. 8p . This conjecture means that the R-string with s = is always unstable 
and a stable R-tube solution, if it exists, must have the following property 

s\ P=0 ^ 0. (3.14) 

In this paper we will assume this conjecture holds and will not consider the constraints from 
the stability of R-string configuration. 

3.3 Rough Estimation for R-tube 

If the domain wall consisting the R-tube is sufficiently thin and resides at p = R, its total 
energy E per the length Az can be estimated as 

1 V(1)(A 2 - R 2 ) + T wall i? + log ( 4V (3.15) 



m 2 27rA2 2 wv ' V(l) 

as has been done in Sec. 12.21 See FigfTl Note that the total energy has divergence terms with 




-^tube ^max 

Figure 7: Potential for a transverse size of R-tube with a rough estimation. 



an IR-cutoff A proportional to A 2 and log A. The former is the energy density V(l) — 1 — A 
of the SUSY breaking vacuum, and the latter is that for the well-known global string tension. 
Here, T wa n is a (dimensionless) tension of the domain wall. This potential has a local minimum 
(maximum) at p = -Rtube(-Rmax) with 

2n 2 n 



V(l) (Vwall - V^all - W 
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if the dimensionless tension T waU is sufficiently large as 

T waU > In. (3.17) 

This is, therefore, a necessary and sufficient condition for existence of the R-tube as long 
as the approximation Eq. ( 13 . 1 5[) is valid. Such configurations that satisfy the inequality 
R > -R max can not avoid to spread out toward the infinite of the space. Note that comparing 
a thickness L wa \\ of the domain wall with -Rtube, if L wa n <C -Rtube holds, the above estimation 
(13.15P works well and we will observe the SUSY vacuum inside the R-tube, namely 



Ip=o 



(3.18) 



Using an approximation discussed in the next subsection, we can show the lower limit of the 
ratio 

^ > n ' > n ' > Vl. f 3 19) 

L wall V(l)T wall L wall 2(1 + e 2 V(l)) 6 ' V ' ; 

Therefore -R tu be can not be very small. If -Rtubc is comparable with Z/ wa u, a configuration of 
R-tube approaches one of the R-string, but s\ p=0 keeps non- vanishing even there in almost all 
the cases as we discussed. 

3.4 Linear Approximation for the Domain Wall 

In order to estimate the transverse size of the R-tube and its stability following the discussion 
in the previous subsection, we need the data {T wa n, L wa u}. We here evaluate them assuming 
that the domain wall in the R-tube can be well approximated by a flat domain wall inter- 
polating the SUSY vacuum at x — — oo and the SUSY breaking vacuum at x = oo. Let us 
consider this configuration in the following. Note that, however, there is an ambiguity for 
definition of T wa u and profile functions for the domain wall since the flat domain wall itself is 
unstable. It is natural to set a relation between the total energy of the system E wa \i and the 
tension T wa u of the domain wall sitting at x — (x) with IR-cutoff A± as 

£waii = / + dx(K + V) = T wall + (A+ - (x»V(l), (3.20) 

which gives a force (pressure) from the SUSY vacuum to the domain wall 

-^£ = V(1)>0. (3.21) 

d {x} 
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Moreover, it is natural to require for the relation, 



dxK 



dxV - (A+ - (x))V(l), 



(3.22) 



is hold near the domain wall solution. When V(l) = holds, the above relation can be derived 
from Derrick's theorem [28]. Then, we define the tension T wa u and a position (x) of the wall 
in terms of only kinetic terms K without using a potential V as 



wall 



dxK, (x) 



dxxK. 



(3.23) 



^wall J\_ 

The equation ( 13 .22 j) is enough to estimate data {T wa n, L wa \\} as the following. Let us 
approximate the profile functions (T, s) for the domain wall by piecewise-linear functions as 

x x 



T 



for < x < L 



wall; 



(3.24) 



-^wall -^wall 

and (T, s) = (1, 0) for x > L wa ii and (T, s) = (0, 1) for x < 0. By inserting this approximation 
to Eq. (13.221) we find that the l.h.s (r.h.s) is proportional to L~l n (L wa n). Note that the tension 
of the domain wall can be expressed as 



wall 



dxK+ dxV- (A + - (z))V(l). 



(3.25) 



Minimizing it in terms of L wa n, we get 

-^wall = 



(A(X) 



c2\2 



(A(X) + e 2 )5(A, e) - (A(A) - C(A) + e»/2)V(l) 



wall 



(A(A) + e») S(A, e«) - ( A(A) - C(A) + - ) V(l) 



(3.26) 



where 



A(X) 



1 

o V(i) 



7T 



7 
60A 

1 



+ 0(A" 



for A > 1 



I V4A - 1 



log(4A - 1) + 0(1) for A 



B(A,e) = - 8- 
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233 1 
l2~4A 



(3.27) 



r 1 1 
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1 xdx arccotA/4A — 1 
2A- 



o V(x) 



V4A- 1 



2 12A 

7T 



+ £>(A~ 2 ) for A > 1 



4V4A - 1 



(9(1) for A 
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Figure 8: Stability of the R-tube with winding number. In the white region, all tube solutions 
n > 1 are unstable. In the purple region, R-tube with n = 1 is stable but others n > 2 are 
unstable. In the light purple region, R-tubes with n = 1,2 is stable but others n > 3 are 
unstable. The region below the red line represents (I3.10p . The SUSY breaking vacuum is 
unstable in the yellow region (see (13.81) ). 

Here we took A_ < and A + > L wa u. Especially we find inequality 



We have used this for deriving the inequality in Eq. (I3.19p . 

Finally, using the result T wa n in the linear approximation, we show the stability condition 
(I3.17P for the R-tube with winding number n in Fig |SJ 

3.5 Numerical approach 

3.5.1 Numerical calculation for stability 

In the previous subsection, exploiting linear approximation, we found the stability condition 
of the R-tube with winding number n. Here, we try to check the parameter dependence 
of the stability by numerical calculation. We adopt a kind of relaxation method to find 
a configuration of the R-tube. See Appendix |X] for details. Since we are interested in a 
parameter region close to borders of the stability of two winding numbers, so we have to treat 
relatively unstable configuration, which require careful analysis. Because of this complexity, 
we focus on a couple of examples for the numerical analysis. 




(3.28) 
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As a first example, we take a parameter e = 1, A = 27/100 where according to the linear 
approximation, winding number n — 1 is stable but n = 2 is unstable (see Fig [8]). Following 
the relaxation method, we take an appropriate initial function and finite relaxation time r, 
then we calculate minimum energy configurations. As we show in Fig and Fig [TU1 energy 
convergences of the configurations have a clear difference in two cases. Here we removed a 
contribution E vev of the vacuum energy density from the total energy E and calculated the 
following dimensionless energy 



and we take the IR-cutoff of the energy as A = 50. The configuration with n = 2 is monoton- 
ically loosing the energy and in a sufficiently late relaxation time r, the energy decreases as 
a linear function, which clearly suggests instability of this configuration. On the other hand, 
as for the configuration with n — 1, the energy seems to converge to a constant value. This 
sharp difference nicely matches with the result of the liner approximation. 

However, it is worth noting that our numerical analysis is done with a finite precision 
which is appropriately chosen by reasonable calculation time. Thus, beyond our calculation 
precision, there may exist an unstable mode which may yield slight energy loss. Thus, as long 
as we use a kind of relaxation method with a fixed initial condition, it may be, in general, 
hard to conclude complete stability of the configuration. However, even if the small instability 
exists, the life-time of R-tube can be longer than the decay time of the R-tube originating from 
an explicit U{1)r breaking effect. In many phenomenological models, the global R-symmetry 
is already broken by adding gravity due to the constant term in superpotential. Thus, at a 
point of the early universe, R-tubes disappear by generating axion domain walls [3Tll39Hi3] . 
Therefore, as long as the stability is long enough compared with its lifetime, we can treat the 
R-tube as a stable solution. 

As a second example, we choose e = 1/50 and A = 6/10. Small e is favorable in model 
building, partially because the longevity of the false vacuum. So from phenomenological 
point of view, it is important to study a configuration of R-tube with small winding number 
in this parameter region. Again, using the relaxation method, we numerically calculate the 
energy convergence of two cases, n = 1 and n = 2. As shown in Fig HU the total energy 
of R-tube with n = 2 converges to a constant value. Thus, within our calculation accuracy, 
the configuration looks stable. Also, the confituration with n = 1 is similarly stable. With 
these numerical results and linear approximation shown in the previous subsection, it may be 
plausible that in small e region R-tubes are relatively stable and the roll-over process does 
not occur. 




A E-K 



A E 



1 



(3.29) 
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Figure 9: Energy against the relaxation time t. Energy of the configuration with winding 
number n = 1 for e = 1 and A = 27/100 converges to a constant value, r is the relaxation 
time. 



Figure 10: Energy against the relaxation time t. Energy of the configuration with winding 
number n = 2 for e = 1 and A = 27/100 decreases linearly at a later relaxation time. 
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Figure 11: Energy against the relaxation time r. Energy convergence of the configuration 
with winding number n = 2 for e = 1/50, A = 6/10. 
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3.5.2 Effective potential for light mode 

As emphasized above, there may exist a very light mode which may cause instability of a 
configuration. Although treatment of such light modes in the relaxation method is not an 
easy task, but we would like to propose a method to uncover the existence of a light mode. 
The most interesting mode is a fluctuation of the size of the R-tube. Generally speaking a 
zero mode (moduli) around a solution is frozen in the relaxation method, and a light mode 
of which dependence in the total energy is quite small, seems to be almost frozen even if 
it exists. To detect such light mode and search the true stable solution, we need to take a 
lot of different initial conditions for the relaxation method. To be concrete, we show initial 
conditions for the fields T, s, % 

/(p) = 1 + tanh( 2 2(P " P0)) tanh(p), (3.30) 

1 - tanh(2(p - p )) 1 + tanh(2(p + p )) 
h V> = 2 2 • 

With various values for po which roughly indicates a transverse size of R-tube, we calculate 
minimum energy configurations with finite relaxation time. For any value of po, the energy 
converges like Figs. [9] and Qj] for those values of n, e and A. However, final configurations can 
have small differences of the energy and the tube-size. To represent the size of the tube, it 
would be useful to introduce the following definition similar to ( 12. 9p . 



o _ J °°dpp 2 K T J °°dpp 2 s'(p) 2 

U T — poo j — r^-j -Ks = f oo , TTyT- {o.oi) 

J dppK T J dpps'{ P y 

Here we defined two sizes of the R-tube, Rt and R s . A reason for introducing two sizes can be 
seen in a discrepancy between the linear approximation in section 13.41 and numerical results 



shown below. As has mentioned in section 12.11 these quantities do not include the cut-off 
dependence and are well-defined candidates for the size of the tube. 

As an example, we take A = 27/100, e = 1. Varying the initial position po, we calculate 
the minimum energy configuration with finite relaxation time. First of all, we show a corre- 
spondence between the initial condition p and the size R s in Fig [121 R s is evaluated with a 
converged configuration. Since it is one-to-one correspondence, varying the initial condition 
Po represents varying the size of the R-tube. 

Now let us show the low energy effective potential for the fluctuation mode of the size. We 
plot the energy (I3.29|) at the relaxation time r = 50, which we will denote as T tu b c = 2ttE(t = 



3 We set Dirichlct condition /(0) = and Neumann condition h'(0) = for T and s for respectively. We 
need to be sensitive for consistency between initial conditions and boundary conditions at p = 0. 
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Figure 12: One to one correspondence between the initial condition po and the size R s for 
e = 1, A = 27/100. 



^tubc 

32.6605 f" 



32.6600 



32.6595 



lo Sio loglTmbc] 

L dr 



-8.0 
-8.5 
-9.0 
-9.5 
-10.0 
-10.5 



Figure 13: Monotonically decreasing potential of R s with A = 27/100, e = l(the left panel). 
T tube = 2tiE(t = 50). Gradients in the energy with respects to r at r = 50 (the right panel). 
Finite gradients indicate instability of conhgurations and noises for small R s imply the limit 
bound of precision of calculations. 



50), with respect to the position of R s in FigJTBl Surprisingly, we observe a monotonically 
decreasing potential in terms of R s for the model with A = 27/100, e = 1 which we explained. 
A large plateau with tiny gradient in FigfTH] is consistent with FigJU] where R s can be almost 
regarded as a massless mode. FigfTS] implies that the R-tube in this case is unstable and will 
expand to the infinity. This clearly suggests that the border line of instability of minimum 
winding R-tube shown in Fig [8] does not match with the numerical analysis above. We find 
that for small R s the two sizes behave differently as shown in FigJT4]and this fact tells us why 
the estimation ( 13.151) with a single size shown in section l3~4l does not work for small R s . It 
would be interesting to study two-scale linear approximation for better approximation. 

Here, we proposed a way to analyze the low energy effective theory for a very light mode 
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Figure 14: Relation between the two transverse sizes Rt,R s with A = 27/100, e = 1. Rt is 
decreasing for small R s whereas Rt for large R s is proportional to R s . Therefore Rt is not a 
good quantity for parametrizing the potential. 



corresponding to the fluctuation of the tube-size and checked the instability of the mode. 
However, this light mode is significantly affected by various corrections such as thermal effects, 
supergravity effects and quantum corrections. We would study these corrections elsewhere. 

4 Bamboo solution: Tube junction 

In the previous section, we showed that an R-string can form a tube-like domain wall with 
winding number and inside of the wall is in the SUSY preserving vacuum. It would be 
wonderful if we could extract any evidences of the existence of the SUSY vacuum from outside 
of the tubqj. Toward this goal, in this section we demonstrate that quantum number in the 
SUSY vacua significantly affects the shape and the tension of the string. Concretely, we 
construct a junction of the R-tube. To the best of our knowledge, this kind of soliton has 
not been known so far. In order to demonstrate an explicit solution, let us again take the 
two-scalar model in Eq. (13. 6ft . As shown in Eq. (13. lip , reflecting the Z2 symmetry of the 
model, there are two different R-tubeJf|. The one has 7 = +1 and the other has 7 = — 1. The 
skin of the R-tube, namely the profile of T field, is independent of the choice of 7, so that one 
can naturally imagine that the two R-tubes with different 7 can be smoothly connected. The 



4 The global R-symmetry is explicitly broken when gravity is coupled to the theory. In this case, when 
the Hubble parameter H becomes the mass scale of the i?-axion, a domain wall interpolating the strings 



is generated and string and walls disappear [3TH39H43] . Thus, one cannot observe a global R-string in the 
present age. However if we replace the R symmetry with another local U{\) symmetry, then a similar tube-like 
solution existing in the present age can be generated. 

5 This is reminiscent of the monopole junction of two cosmic strings studied in [44] . It would be interesting 
to study further our junction in light of this similarity. 
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Figure 15: The bamboo soliton (A slice at the center): Junction of the R-tubes with 7 = +1 
and 7 = — 1. In the left panel, the red region is where s > 1/2 while the blue is the region 
where s < —1/2, and the yellow stands for the region where 7/20 < \T\ < 12/20. The right 
panel shows the potential density isosurface with which one can clearly recognize the domain 
wall inside the tube. 

junction of the R-tubes is a domain wall which interpolates two different SUSY vacua inside 
the R-tube. We call this junction the R-bamboo. 

A numerical solution for e = 1, A = 27/100 is shown in Figs. [15] and [161 Far away from 
the domain wall along the tube, the solution asymptotically goes to the R-tube solution. At 
the junction, the transverse size of the tube becomes smaller since the domain wall pulls the 
tube toward its inside, see Fig. [T5l 

This R-bamboo solution may be created when two R-tubes collide. If the two tubes are 
different kind, the domain wall must be created at the junction of the two tubes. At the 
same time, the anti R-bamboo may be created. This is very similar phenomenon to monopole 
and anti-monopole creation associated with the non-Abelian string reconnection. This is 
interesting issue but is beyond the scope of this paper, so we leave it as a future work. 

Finally, it is worthy to note that stability of bamboo configuration is not guaranteed by 
our numerical approach. As emphasized in the previous section, our analysis is a kind of 
relaxation method. Actually, if there is no domain wall inside the R-tube, the configuration 
seems to have a small instability shown in Fig [131 However, in the bamboo configuration, 
the existence of the domain wall increases the stability because the energy of the domain 
wall is proportional to R 2 . Thus, if the number of the domain walls is large, the bamboo 
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Figure 16: 3D plot of the Bamboo soliton for e = 1, A = 27/100. 

configuration would be stabilized enough. In principle, the same argument done in Fig [13] 
would be applicable, however numerical calculation becomes quite involved because one needs 
the two-dimensional relaxation method. 

5 Conclusion 

In this paper, we have demonstrated a fascinating role of a cosmic R-string/R-tube by using 
two toy models with spontaneous R-symmetry breaking. The first example shown in section 
12. II is a single field model. The model can be regarded as a toy model of one of O'Raifeartaigh 
models studied in [12] in which pseudo-moduli space is stable everywhere. In this model, 
string-like defect generated by the Kibble-Zurek mechanism is stable and very close to the 
known global strings. Winding number dependence of the size of the string is linear in n. On 
the other hand, in the two-field model shown in section [31 a string-like object is a tube-like 
domain wall interpolating a false vacuum and a true SUSY vacuum. One specific feature of 
two-field model is existence of a tachyonic direction at R-symmetry restoring point X = 0. 
Because of this, the core of the string is not stable under fluctuation toward the tachyonic 
direction. Thus, inside the core in which R-symmetry is restored, can be filled out by the 
true vacuum through the tachyonic direction. Naively, one may think that such an R-tube 
is unstable. However, as we shown in the main text numerically, there exist metastable R- 
tube solutions for certain parameters. An interesting property of such R-tube is the winding 
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number dependence of the size. By the linear approximation, we estimated the dependence 
and found that its dependence is n 2 rather than n. Therefore, there is a tendency that tube 
configuration with larger winding number is more unstable. Numerically, we checked the 
higher winding instability at a sample parameter e = 1 and A = 27/100: We showed that 
the configuration with winding number n = 2 is unstable and the total energy decreases 
monotonically. 

If an unstable R-tube is created by the Kibble-Zurek mechanism, it rapidly expands and 
our universe will be completely filled by the true vacuum. This process gives constraints 
for models building. However, it is worthy to emphasize that such roll-over process can be 
protected by D-term contribution or thermal potential. As is demonstrated in [2TIB5| B6]. 
when a D-term contribution cannot be negligible, it can lift the tachyonic direction and 
stabilize the pseudo-moduli space. In such models, the roll-over process does not occur. Also, 
when the amplitude of (tachyonic) messenger mass at the origin is sufficiently smaller than 
that of R-symmetry breaking field, the vacuum selection is successfully realized by exploiting 
the thermal potential or Hubble induced mass. If such thermal potential keeps lifting the 
tachyonic direction until the time of the R-string decay [3T |[39l4^3| . the roll-over process can 
be successfully circumvented. However, even if such an early stage scenario is assumed, there 
are sever cosmological constraints on R-axion density as studied in [3Tj . 

Finally we comment on our numerical analysis done in section[3l Since we adopted a kind of 
relaxation method with a fixed initial condition to find energetically minimum configurations, 
it is not easy to conclude the stability of the configuration sharply. There may exist a very 
light mode which does not change the energy significantly. To uncover the existence of the 
very light mode, we proposed a method for studying the low energy effective potential for such 
light mode. By changing the initial conditions and observing converged configurations, we 
can estimate the effective potential. Using this method, at the parameter e = 1 A = 27/100, 
we find emergence of a light model and find a large plateau in the effective potential. It would 
be important to apply the method to a wide range of parameter space and study borders of 
instabilities of configurations with various winding numbers numerically. This is beyond the 
our scope, so we will leave it as a future work. 
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A Relaxation method 

In this section we review a kind of relaxation methods applied to constructing numerical 
solutions in this paper. Let us consider the following general Lagrangian for scalars Q in 
~R d+1 with coordinates = {t,x 1 }, 

C=^g a ^ a d^ -V(cf>), (A.l) 

where g a p is the metric for the field space. Our goal is to find an numerical static solution 
4> a (t, x l ) = 4>s Q \{x l ) for this system. The shooting method is a good strategy for a system with 
a single scalar fields and a single spatial coordinate x , but only in that case. For systems 
with multi fields in higher dimensions, the shooting method does not work very well and we 
need the relaxation method explained bellow. Let us introduce a 'relaxation time' r instead 
of the real time t and suppose that 4> a depend on r as (p a (r, x l ). Then r dependence of cf> a is 
defined by, 



5C 8C 



with Neumann condition at the boundary of the region E C 



= 0. (A.3) 



The added term in the r.h.s of Eq. (lA.2|) i the so-called friction term. Actually, due to this term 
we can show that the ordinary total energy E with integral region E is no longer constant 
but a monotonic decreasing function E = E(r) of r as 

dE(r) 



dr J s dr dr 

If we observe the energy converges, we find a solution 0" ol as 



d d xg a p < 0. (A.4) 



lim EM = E sol & lim X' X ] =0 <S> lim <j> a {T,j) = Ci(^)- (A.5) 

t— >oo t—>oo OT t— >oo 
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If the energy has the global minimum, therefore, we can obtain, at least, one static solution 
0" ol by using the relaxation method. If there exists multiple local minima of the energy, an 
initial condition of <p a chooses one of them. 

In actual numerical computation we need a cutoff of the relaxation time at r = Tfi n 
and we regard d> a ( solution. A relation between Tfi n and precision of a solution 

<j) a = a (rfj n , x l ) can be discussed in the following. If deviations of (f) a from a solution </>" ol are 
efficiently small, they can be expanded as 

oo 
n=l 

with profile functions f n {% % ) for massive modes of mass m n around the configuration 0" ol . 
Eq. (1A.2p gives development of coefficients a n (r) as 

a n (r) = <e-"^, (A.7) 

and behavior of the energy is controlled by the lightest mass mi as 

E(t) « E sol + Ae~ 2m " T . (A.8) 

Here a constant A G M>o depends on an initial condition we took and is assumed to be the 
same order as -E so i. To get precision 10~ p , therefore, we have to take a time r^ n for this 
relaxation method as 

r fin > In 10, (A.9) 

where mi can be roughly guessed by a typical mass scale of the system. With large r, we 
sometimes observe a random behavior of E{r) which is a signal of the limit bound of machine 
precision. See FigJT7]for an example. 

If a configuration of <p a is accidentally near to a saddle point, we also observe an exponential 
decay of E(t), but after that it collapses like a waterfall as 

E(r) « £ saddle + Ae~ 2m2r - Ae 2 ^ T (A.10) 

with a tachyonic mass m 2 = — \m 2 \. Therefore an exponential behavior of the energy do not 
always guarantee that a stable solution is obtained. Taking multi initial conditions of 4> a can 
avoid this technical error. 
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Figure 17: Typical behavior of an energy in the relaxation method. Here we get accuracy of 
10 -7 ' 5 for a solution. Calculations after r ~ 32 turns out to be meaningless. 
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